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We construct initial data corresponding to a single perturbed Kerr black hole in vacuum. These 
data are defined on specific hyperboloidal slices on which the mean extrinsic curvature K asymptot- 
ically approaches a constant at future null infinity . More precisely, we require that K obeys the 
Taylor expansion K = Ko + 0{a'^) where Ko is a constant and a describes a compactified spatial 
coordinate such that is represented by a = 0. We excise the singular interior of the black hole 
and assume a marginally outer trapped surface as inner boundary of the computational domain. 
The momentum and Hamiltonian constraints are solved by means of pseudo-spectral methods. We 
find exponential rates of convergence of our numerical solutions, and plan dynamical evolutions in 
a future project. 
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I. INTRODUCTION 

One particular objective of numerical relativity is the 
study of gravitational waves emitted by the interaction 
of astrophysical objects. The corresponding wave forms 
permit the construction of templates for gravitational 
wave detections. Most of these calculations have been 
carried out on spatially truncated numerical domains. 
Despite their remarkable successes in recent years, the ac- 
curacy obtained was limited due to effects coming from 
artificial outer boundary conditions. Moreover, as the 
outgoing gravitational radiation is well-defined only at 
future null infinity J^^ , the wave extraction at finite dis- 
tances causes additional inaccuracies. 

An elegant alternative approach, which avoids these 
undesired effects, is the inclusion of as outer bound- 
ary of the computational domain, i.e., the study of the 
spacetime on compactified hyperboloidal slices (see [I 
for a review) . This description permits the detailed anal- 
ysis of the outgoing radiation through the Bondi news 
function [21 E]. 

Note that, in this approach, the physical metric be- 
comes singular at J^^ , which gives rise to a decomposi- 
tion into a conformal factor, absorbing the singular be- 
havior, and a regular conformal metric. Friedrich [?] has 
formulated a hyperbolic system for the conformal metric 
on hyperbolidal slices, which is equivalent to the Einstein 
equations. Unlike the Einstein equations, this system is 
manifestly regular at J^+. The first numerical results of 
this system (in the context of weak relativistic fields were 
obtained) by Hiibner [5]. 

Moncrief and Rinne [B] showed that the apparently 
singular boundary terms appearing in the hyperboloidal 
concept can be explicitly evaluated at J^"'" in terms of 
conformally regular geometric data. As a result of this 
analysis, they obtained a rather rigidly constrained and 
gauge fixed formulation of the field equations. The first 
stable dynamical numerical evolution based on this for- 
mulation was published by Rinne 7 . The corresponding 
code provides long-term stable and convergent evolutions 



of axisymmetrically perturbed Schwarzschild spacetimes. 

The construction of initial data as a starting point for a 
subsequent dynamical evolution of the space time geom- 
etry requires the solution of coupled elliptic constraint 
equations, specifically the Hamiltonian and momentum 
constraints. In this paper we consider data which are 
defined on spatially compactified hyperboloidal slices ex- 
tending up to Previous calculations of such data 
have been carried out in the context of constant mean 
curvature ( "CMC-condition": K = constant throughout 
the slice) . Important contributions were made by Frauen- 
diener [S] and Buchman et al. While Frauendiener 
looked at slices with isotropic extrinsic curvature (that is, 
the extrinsic curvature is assumed to be proportional to 
the conformal 3-metric in the slice) , Buchman et al. com- 
puted data in the context of conformal flatness, taking 
a generalization of the Bowen-York curvature |10j as the 
solution to the corresponding momentum constraint. In 
these cases, only the Hamiltonian constraint needs to be 
solved, and it is degenerate at J^"*" by virtue of the sin- 
gular physical metric. For isotropic extrinsic curvature, 
the Hamiltonian constraint reduces to a form which is 
sometimes called the Yamahe equation. This equation 
has been studied extensively and solved numerically with 
pseudo-spectral methods by Frauendiener [S]. Anders- 
son proved conditions under which the Yamabe equation 
has regular solutions [TT]. These criteria were general- 
ized by Andersson and Chrusciel [T5J [T3] to the situation 
in which arbitrary slices obeying the CMC-condition are 
considered. As an example, these mathematical results 
apply to the slices used by Buchman et al. [3]. 

In this paper, we construct axially symmetric, confor- 
mally non-flat initial data corresponding to single per- 
turbed Kerr black holes in vacuum on hyperboloidal 
slices. The corresponding constraint equations are con- 
sidered in the conformal transverse traceless decomposi- 
tion '141 115j. In this formulation, the momentum con- 
straint contains a term behaving like a^^V^K , where a 
describes a compactified spatial coordinate such that ^+ 
is represented by cr = 0. Now we can regularize the mo- 
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mentum constraint by requiring that the mean curvature 
K obeys the Taylor expansion K = KQ+0{cr^) where Kq 
is a constant. Henceforth, we call such sUces "ACMC"- 
slices, as they fulfil the condition of asymptotically con- 
stant mean curvature. It turns out that the degenerate 
behavior of the corresponding Hamiltonian constraint is 
asymptotically the same as that of the Yamabe equation. 
For this reason, we call this behavior in the following 
"Yamabe-like." 

As a starting point for our computations we take the 
Kerr solution in Kerr coordinates which are horizon- 
penetrating. Using the "height- function" -technique from 
Zenginoglu '16j, we apply a particular coordinate trans- 
formation in order to obtain compactified hyperboloidal 
ACMC-slices. This is to say, the coordinate transfor- 
mation is chosen such that, in these compactified slices, 
the extrinsic mean curvature K is asymptotically con- 
stant, K = Kq -\- 0{a'^). Now we consider the regular 
momentum constraint as well as the Yamabe-like Hamil- 
tonian constraint that arises in the conformal transverse 
traceless approach. The corresponding computational 
domain is bounded by (outer boundary) and by a 
marginally outer-trapped surface which serves as inner 
boundary. In other words, we excise the black hole in- 
terior and require that at the excision surface, the ex- 
pansion of outgoing null rays vanishes. Of course, for 
the Kerr solution we started with, the constraints are 
identically satisfied. It is, however, easy to modify the 
set-up slightly and to compute initial data which cor- 
respond to perturbed Kerr black holes. To be more spe- 
cific, we keep the coefficients, appearing in the constraint 
equations, as well as the boundary values of the solution 
to the momentum constraint as they are given by the 
Kerr solution in our ACMC-slices, but modify the coor- 
dinate location of the marginally outer-trapped surface. 
We then solve the constraints in the new computational 
domain given by the modified inner boundary and ob- 
tain perturbed Kerr black hole data. The deviation from 
the (unperturbed) Kerr solution is determined by a sub- 
sequent analysis. This analysis consists of the compu- 
tation of mass and angular momentum multipoles that 
characterize the marginally outer-trapped surface. We 
find distinct differences with respect to the correspond- 
ing Kerr expressions which means that our numerical so- 
lutions cannot describe an unperturbed Kerr black hole. 

The computation of these initial data has been real- 
ized by utilizing a single domain pseudo-spectral method 
which provides us with a rapid exponential convergence 
rate of the numerical solutions. This exponential fall-off 
is a strong indication that for the ACMC-slices the same 
regularity criteria hold as for strict CMC-slices. 

The paper is organized as follows. In section |ll] we 
construct hyperboloidal ACMC-slices of the Kerr solu- 
tion. Next, in section |TTT] we discuss the constraint equa- 
tions within the conformal transverse traceless decompo- 
sition and elaborate on regularity criteria valid at . 
The inner boundary is the subject of section |IV| while 
we present the conditioning as well as the description 



of the numerical treatment is presented in section [V] In 
section |VI| we provide exemplary results corresponding 
to strongly perturbed Kerr data as well as to weakly 
perturbed data with large rotation parameter. We con- 
clude in section |VII| with a discussion and a perspec- 
tive of future work. The appendix contains the explicit 
form of the auxiliary function A used in the coordi- 
nate transformation to obtain ACMC-slices (see section 



VIII A I, a comparison between CMC- and ACMC-slices 



in the Schwarzschild case (see VIII B ) and, finally, section 
|VIII C| in which the multipole analysis of the constructed 
data is described. 

Note that we use the following conventions. Greek 
indices run from to 3 where denotes the time com- 
ponent. Small latin indices are from 1 to 3, covering 
spatial indices. Conformal (unphysical) objects will be 
decorated with a tilde. We use units in which the speed 
of light as well as Newton's constant of gravitation are 
unity. 



II. THE KERR SOLUTION IN ACMC SLICES 

As mentioned above, the construction of initial data 
pursued in this paper is based on the conformal trans- 
verse traceless decomposition of the constraints on com- 
pactified hyperboloidal slices. In order to regularize the 
momentum constraint in this formulation, we consider 
ACMC-slices on which the mean curvature K obeys the 
Taylor expansion K — Kq 0{a^). In the following, 
we identify such slices within the Kerr solution, which 
then will serve as a starting point for the generation of 
perturbed Kerr data. 

We begin with the horizon-penetrating line element in 
Kerr coordinates {V, r, 9, (p), 
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Here we have 
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where M and J = aM are, respectively, the mass and 
angular momentum of the black hole. Note that in the 
following we use the dimensionless angular momentum 
parameter 
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The event h orizon of the black hole is located at r = 
Th — M(l -I- v^l — p) which runs from rh = 2M in the 
Schwarzschild case {j — 0) to = M for the extreme 
Kerr black hole = 1). 
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We introduce hyperboloidal slices through a coordi- 
nate transformation that is motivated by an asymptotic 
integration of outgoing radial null rays, defined by 

ds^ = 0, de^ 0, dLp = 0, dVdr > 0. 

The corresponding equation 

2dr 



dV 



1 



2Mr 

r,2 



1 



2M 

r 



0{r- 



dr 



(3) 



is solved asymptotically by 



constant 



f — 



log 1 



From this form, we obtain a coordinate transformation 
which leads us from Kerr coordinates to general axi- 
ally symmetric, compactified hyperboloidal slices by the 
ansat^ 

r= (4) 



V = AM 



log((T) + yl(o-, cos( 



(5) 



In the new coordinates (r, cr, 0, (p), the compactified hy- 
perboloidal slices are described by r = constant. The 
coordinate location of future null infinity is given by 
a — Q while the event horizon is placed at 



cr = (Th = 



1 



2 



(6) 



The coefficients g^^ of the physical metric can be written 
in the form 



(7) 



Here, g^^ denotes a conformal metric, which is regular 
everywhere, and is a conformal factor that vanishes 



linearly at , i.e. Q,\ 



0, d„n\jf+ ^ 0. For the 



function A appearing in (|5|, an arbitrary regular expres- 
sion in terms of a and 

IJL := cos 6* e [-1, 1] 

can be choserj^ as long as one ensures that the result- 
ing hyperboloidal slices r = constant are spacelike every- 
where away from . This construction of hyperboloidal 
slices is equivalent to the height function technique |16| . 
and the corresponding height function h is indicated in 



In the following we make use of the freedom to choose 
the regular function A in order to generate slices that 
satisfy the ACMC-condition 



K = Ko + 0{<7^), Ko^ constant . 
Here K is the mean curvature 

K = /^^,.g''^ 
i.e., the trace of the extrinsic curvature K^^, 



(8) 



(9) 



(10) 



In this definition, £„ is the Lie derivative along the future 
pointing unit vector normal to the slices r = constant. 
We follow the sign convention employed by Wald |T7], 
which was also used in [9 . In this convention (opposite 
to the widely adopted one by Misner, Thorne & Wheeler 
|18p. the mean curvature K assumes positive values at 

Now, for a given function A(cr, /i), the Taylor expansion 
of K with respect to a at {a — 0), considered in a 
slice T = constant, has the following structure: 



Kq + Ki(T^ + K2cr^ + Ksa^ + ©(a"*) 



(11) 



with rather longish expressions for the coefficients Kj — 
Kj{ii). Putting Kq = constant and Ki = K2 — K3 = 0, 
we obtain asymptotic conditions for A(a,fi), 

= ao 



A 

d,A 

dlA 
dlA 



a=0 



<T = 



ai 



16 

o2 



-2(4 + 6ai -(-3a?) 



(12) 

(13) 

(14) 
(15) 



+j' {l+-a'{l + a^) 
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with flo and ai being arbitrary, unconstrained constants. 

As mentioned above, we need to ensure that the corre- 
sponding slices r = constant are spacelike, which means 
that the conformal lapse 

1 

a = ria 



is well defined everywhere. Now, for numerical reasons, 
it would be desirable to have a "very" smooth function 
(5, meaning that a spectral expansion with respect to 
Chcbyshcv polynomials shows a rapid exponential fall-off 

of the corresponding coefficients c^"-* Indeed we ffnd 



^ Note that through this ansatz we retain the horizon penetrating 

nature of the coordinates. ^ As a representative, we consider a{a-, fj, = 0), i.e. the conformal 

^ The requirement of axisymmetry imphes that A has to depend lapse as a function of cr £ [0,o"h] along the equatorial plane 9 = 

on cosd, which we denote henceforth by /i. ■7t/2. 
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that a rapid convergence rate of these coefficients is a 
good indication of high numerical accuracy of the results 
obtained from the subsequent pseudo-spectral solution of 
the constraints. We thus aim at a choice for the function 



A such that on the one hand the conditions ( 12 1-( 15 1 are 
satisfied and on the other hand the jc^"' | fall off rapidly. 

We have found that these requirements are met by an 
appropriate diagonal Pade expansion. In this approach, 
the coefficients bk = bk{ij) and Ck — Ck{lJ.) of the rational 
function 



E ha'' 

71 

1 + E ck<j' 

k=l 



(16) 



are determined by the first 2n derivatives of A at cr = 0, 
i.e. through the conditions given in (12l-(15). Since oq 



represents an irrelevant time shift, we can set it to zero, 
which implies 5o = 0. In order to satisfy the remaining 



three conditions ( 13 1-( 15 ), we need at least a second order 
Pade expansion, n = 2, which provides us with the free- 
dom to choose a parameter, say 04, describing the fourth 

derivative d^A .In this way, our Pade approximation 

(T = 

contains two undetermined parameters, ai and 04, and 
we can use this freedom in order to obtain rapid fall-off 
rates of the corresponding spectral coefficients c^"-*. A 
satisfactory choice is given by 

113 9 , ., 3 



ai 



04 



100 500 



4! 



161 

To" 



5 .2 



(17) 



(18) 



which describes the two constants parametrically in 
terms of the dimensionless angular momentum j. The 
explicit form of the resulting function 



A{a,^l) = ^2{(T,^l) 



(19) 



as well as a discussion regarding the derivation of the 
choice ( |T7|l8l ) can be found in Appendix IVIII AI Note 
that the corresponding mean curvature K assumes the 
required ACMC-structure, 



X(a,/i) 



M^16{l + a,)-f 



(20) 



III. THE CONSTRAINTS IN THE 
CONFORMAL TRANSVERSE TRACELESS 
DECOMPOSITION ON COMPACTIFIED 
HYPERBOLOIDAL ACMC-SLICES 



with a being the lapse function, the shift vector, and 
7ij the induced three-metric of the spatial slice. The 
constraint equations valid in that slice are the momentum 
constraints 



V,(if*^' - 7*^' if) = 
and the Hamiltonian constraint 



(22) 



R+ - K,jK'^ = 0. (23) 

Here, Vi is the covariant derivative associated with the 
induced metric 7^, and R is the associated curvature 
scalar, while the Kij denote the spatial components of 
the extrinsic curvature. 

There are several formulations for decomposing the 
constraint equations. Here we use the conformal trans- 
verse traceless decomposition (CTT) |14l[TBj . in which we 
write 



1 



(24) 



o 

with the conformal three-metric defined as 



7,, = fj- = f]-v^ 



and the traceless part A*^ , 

7.,i'-'-0 = fJi,,. 
We further decompose A^^ as 

where 



(25) 



o 

and Vi denotes the covariant derivative with respect to 

lij- 

In the CTT approach, the symmetric tracefree tensor 
M'-', the mean curvature K, and the conformal metric 
7y- are considered as freely specifiable data. Given these, 
the momentum constraint becomes 



(26) 



with 



A. Constraint equations 



In the standard 3-1-1 decomposition of spacetime, the 
line element is written as 



To rewrite the Hamiltonian constraint, we use the de- 
composition (see, e.g., |T]) 



^dt^ + {dx' + 13'dt) (dx^ + put) . (21) 



(27) 
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where w is a prescribed spatial functiorj^ which describes 
by = 0. Note that a; has similar properties to 

ri, i.e., it is positive away from and possesses non- 
vanishing gradient at The auxiliary potential (j) is 
chosen to be positive everywhere. 

The Hamiltonian constraint then turns into 



(28) 



Here R denotes the curvature scalar associated with 
Note that for isotropic extrinsic curvature, Aij = 0, the 



last term vanishes, and (28) turns into the Yamabe equa- 

Finally, with the decomposition of the conformal factor 
the momentum constraints read as 



(271 



(29) 



In sum, the determination of valid initial data on com- 
pactified hyperboloidal slices in the CTT approach con- 
sists of the solution of the conformal constraint system 
(28 29 1 for the potentials (p and . As mentioned above, 
M*^ , K,^ij as well as w are prescribed quantities^^and A^^ 
is given in terms of and M*^ by equation ( 25 ) . 



B. Regularity of the constraint equations 



The momentum constraint ( 29 1 
■ ^^^V^K which is regular up to 



contains a term 
only when the 



mean curvature satisfies the ACMC-condition 
K — Kq + K4Uj^, Kq — constant. 



(30) 



If the global CMC condition is imposed, that is K4 = 0, 
then the constraints decouple, i.e., at first, the momen- 
tum constraints can be solved for and then the Hamil- 
tonian constraint for (f>. Although we lose this separa- 
bility for K4 ^ 0, we retain the structure of the con- 
straint equations at .y^. To be more specific, the mo- 
mentum constraint is a completely regular elliptic system 
of differential equations, requiring the prescription of in- 
ner and outer boundary conditions for a unique solution. 
At the same time, the Hamiltonian constraint assumes a 
Yamabe-like structure, i.e. it looks asymptotically up to 
order like the Yamabe equation. As a consequence, 
the first three terms 4'o,4>i,4>2 of a Taylor expansion of 
(j), taken in normal direction with respect to , can be 
found explicitly. 



t>2UJ 



higher order terms, (31) 



* In this paper, we choose uj to always agree with the compactified 
radial coordinate cr. 

^ Isotropic extrinsic curvature implies K = constant, i.e., the cor- 
responding slices satisfy the CMC-condition. 



e.g., for 00 we obtain: 



(32) 



All terms higher than quadratic order cannot be deter- 
mined a priori, but follow from the global solution of the 
constraint system if boundary conditions are imposed. 

Andersson and Chrusciel investigated regularity cri- 
teria valid for solutions to the Hamiltonian constraint 
on hyperboloidal slices satisfying the CMC-condition 
[TT| [T^ . They found that cj) (and hence H.) is analytic 
in a vicinity of provided that at the conformal 
extrinsic curvature 



f^ij — 2^^sqij (33) 
is proportional to the induced conformal metric 

qi] = li] - SiSj, (34) 
that is, that there is a scalar A defined on J^'^ such that 



(Rij - Xqij) \y+ = 0. 



(35) 



Here, P is the conformal unit vector normal to ^+ in 



the hyperboloidal slice. The condition (35 1 is equivalent 
to requiring that be shear-free. 

In this work we assume that the shear-freeness of 



( 35 1 ensures likewise the regularity of solutions to the 
constraint system (28 |29) which arise in the context of 
the more general ACMC-condition. As we use for the 
conformal metric 'jij the Kerr expression in ACMC-slices 
(constructed in section |ll|, we can explicitly check the 
validity of (351. We find A = 0, i.e., the condition is 



fulfilled in a trivial manner. 

Note that for CMC slice in general, the Taylor expan- 
sion of the auxiliary potential (j), taken in the normal 
direction with respect to J^'^ , breaks down at 3rd order 
in w, i.e. , terms oc u!^log{uj) appear. These logarithmic 



terms vanish if the shear-freeness of [i.e., Eq. (35 1] 



is imposed. Bardeen et al. [H] discuss the fact that this 
condition is equivalent to the absence of outgoing radia- 
tion at . Here we are interested in completely regular 
initial data which do not possess these logarithmic terms. 
Following the arguments presented in [12], we conclude 
that the data constructed below do not carry outgoing 
radiation. 



IV. MARGINALLY OUTER-TRAPPED 
SURFACES AS INNER BOUNDARIES 

In this paper, we compute initial data corresponding 
to a perturbed Kerr black hole in vacuum. In order to 
describe the black hole, we excise its singular interior and 
require the excision boundary H to be a marginally outer- 
trapped surface [20J. We thus obtain an inner boundary 
of our computational domain. 
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For a marginally outer-trapped surface, the expansion 
scalar of the 2-manifold embedded into the 3-slices van- 
ishes for outgoing null rays. The corresponding condition 
reads 



V„ 



K - s'^s^'Kat = 0, 



(36) 



where s" is the outward pointing unit vector normal to 
the surface H in the 3-slice. 

Note that in ACMC-slices of the Kerr solution (see 
section |ll| , this condition is satisfied at the location of 
the event horizon r — r\^, corresponding to cr = Ch = 
2/(1 + ^/l-f)■ Further more, this surface is the outer- 
most marginally trapped surface, i.e., the apparent hori- 
zon [21]. 



Writing (36 1 in terms of conformal quantities, we get 
an inner boundary condition valid at T-L, containing (f) as 
well as its first derivatives, 



4 \ w / 6a; 



S'^tAab^Q. (37) 



Here we used the equations (|7 24 and 27 1 and introduced 
the conformal normal vector s* = fi^^s*. 

For the construction of perturbed Kerr data below, we 
require this condition at a coordinate location a-u that 
differs from the Kerr value Ch [see Eq. ([6])]. At the same 
time, we keep all coefficients in the constraint equations 
as well as the boundary conditions of the vector in the 
form that results from the Kerr solution in ACMC slices. 
If the perturbation is chosen to be sufficiently small, we 
ensure that the corresponding marginally outer-trapped 
surface will be outermost [22^.23^, i.e., the apparent hori- 
zoijEl 



where the iJ, are defined as: 



Hi 

H2 

H3 
Hi 
H5 
He 

Hj 
Hs 



16^2712 

-8aV''f2,-8af 1 
-^2^ + 4^^"?!, + 67" 

2 , 

3 



(39) 



By virtue of the decomposition ( 25 ) we have written here 



A^^ and Aij as 



ya 



m • 



where £, I, m are given by 



and 



a 

= 

a 

m'^ = 

Uja 



jkif.3 



M,,. 



The momentum constraints read 



(40) 



(41) 



(42) 



V. NUMERICAL SOLUTION OF THE 
CONSTRAINT EQUATIONS 

In this section we explicitly describe the steps in the 
set-up of the constraints and their numerical solution. 
Using the coordinates {x^ , , x'^ , x^) = (r, cr, (/?) and 
putting uj — a, we write the Hamiltonian constraint as 

+ Hj(j)^ + Hscj)"' = 0, (38) 



^ Even though stronger perturbations may lead to solutions for 
which the apparent horizon could be located further out 1241 . 
supplementary investigations of the computed solutions strongly 
indicate that, in these cases, we are also dealing with the appar- 
ent horizon (see appendix | VIII C| l . 



a"%^':,r + b'^pVP^ + t\V^ + 0* + eV^ = 0, 
where 



irq 



1 



(43) 



r^f^^-7^'=f«,5;+7^'fjp+7^'=np 



2 

^3 



kp) 



+1'" (f I 



F!- ,r 



^Pd 3^ Ip ' jl kp) 

^PJ Ip jl kp) 



(44) 



3a3 



In order to derive particular numerical values for this 
set of conformal coefficients, we first need to specify the 
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conformal factor r^Korr that we take for the conformal 
decomposition of the unperturbed Kerr metric in ACMC 
slices. We put 



Kerr 



4M' 



a to ^Kcrr = 2 

Fraktur letters 



M. Then 
well 



leading via (27) with ui 

all coefficients denoted m I'raktur letters as well as 
{Hi, . . . , H^} follow explicitly from the Kerr solution in 
ACMC slices. In particular, the quantities e* are regular 
at J^^ by virtue of the ACMC condition [cf. equation 
(20 1]. The coefficients m*^ are obtained by requiring that 



the unperturbed Kerr solution satisfie s th e momentum 
constraints through Vj^^rr = 0- Using (25 1, we conclude 
that 

Now, for solving the constraints uniquely, we still need 
to impose boundary conditions. As the momentum con- 
straints are a completely regular, non-degenerate system 
of elliptic equations, we may simply choose trivial Dirich- 
let boundary conditions for the potentials V^, 



V = = v 



(45) 



Note that other choices are possible. Indeed any condi- 
tion that is compatible with l^err = in the unperturbed 
case will do. In contrast, the potential is not subject to 
freely specifiable boundary conditions. On the one hand, 
the Yamabe-like Hamiltonian constraint provides us with 
a regularity condition valid at J^^ , 



^97 



K- 



(46) 



which is equivalent to (32 1. On the other hand, we re- 
quire our inner boundary H to be a marginally outer 
trapped surface, which means that (37 1 yields the bound- 
ary condition 



where the t)i are defined through 
1 







n 



(47) 



f)2 = 



f)3 



K 
6a 



(48) 



As already mentioned, we impose this condition at cr = 
a-fi ^ (Jh [see Eq. ([6|], in order to obtain perturbed Kerr 
data. This is the only perturbation we perform, though 
various other choices are possible. We solve the cou- 



pled eUiptic system (38 



ary conditions (45 46 47 



pseudo-spectral Gauss-Lo 



43 1 together with the bound- 
by means of a single-domain 
jatto method. For prescribed 



TABLE I. Sample data for the construction of the initial data. 



Calculation No. 


j 


CTh 




1 


0.5 


1.0718 


0.5 


2 


0.9 


1.39286 


1.4 


3 


0.999 


1.91441 


1.9145 



numerical resolutions Nu 
are given through 



and Nu, the collocation points 



- ^?^sin2 ] (j = 0,...,iV,) 



^fc = -cos 



Trfc 
A^ 



(49) 



(fc = 0,...,iV^), 



which are the abscissa values of the local extrema of the 
Chebyshev polynomials Tiq^{2a / a-u — 1) and T^^^fi) re- 
spectiveljj^ As described in detail e.g., in [25J, the values 
{4'jk^V^J,, V^tj, VjTj) of the unknown potentials at these col- 
location points are combined to form a vector f(^'''^f). 
From any such vector, the Chebyshev coefficients of the 
potentials as well as those of their first and second deriva- 
tives can be computed. The combination of constraints 
and boundary conditions, evaluated at the points (49 1, 



yields a non-linear system of algebraic equations of or- 
der 4 X {N^ -t- 1) X {N^ + 1) for the entries of the vector 
^{N„,N^) ^ We solve this system by means of a Newton- 
Raphson scheme, which uses the iterative "bi-conjugate 
gradient stabilized method" [25] for inverting the Jaco- 
bian. This method works well only if an appropriate pre- 
conditioner is supplied. Here we utilize a finite difference 
representation of the Jacobian and invert it with the help 
of a band diagonal matrix decomposition algorithm (see, 
e-g-, [57] and references therein). An initial guess for the 
Newton-Raphson scheme is taken from the unperturbed 



solution, {(/)Korr = 2V M, Vj^orr = 0} for which CT-H = (Th- 

We depart gradually from this solution by modifying the 
location au of the inner boundary in several steps. Any 
solution computed serves as the Newton-Raphson guess 
for the subsequent solution. In this way we can obtain 
large perturbations, as will be depicted in the next sec- 
tion. 



VI. RESULTS 

For any prescribed rotation parameter j with 
5: bl < li the numerical scheme permits the construc- 
tion of initial data describing perturbations of the corre- 
sponding Kerr black hole. The deviation is tuned through 



^ As our perturbations lead to equatorially symmetric solutions, 
the restricted interval £ [0, 1] would suffice. However, we have 
chosen to keep the full interval to allow for non-equatorially sym- 
metric perturbations that might become relevant in the future. 
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the prescription of the location a-u of the inner boundary. 
Values of sample data are given in Table [ij 

We observe that the sign of the deviation from the 
unperturbed solution depends on the sign of (cr^ — ah), 
i.e., the location of the inner boundary with respect to 
the unperturbed solution. 

In all calculations we find a rapid exponential decay of 
the Chebyshev coefficients of our numerical solution up 
to machine precision. Moreover, we see an exponential 
convergence rate with respect to the numerical resolution 
(ncr,n^) = (iVcr + l,iV^ + 1). In particular, we monitor 
the quantities 



spectral convergence rate 



=M-i/2sup|</>„ 



"60,301 



^60,30 I 



(50) 
(51) 



which describe the maximal deviation from the numeri- 
cal reference solution with {ncr,n^) = (60,30). For ap- 
proximating the supremum, we take coordinate values 
cr G [0,(t-h], /i € [—1,1] located on an equidistant grid 
of 500 times 500 points. For the examples listed in the 
table, we provide in figure [T] logarithmic plots displaying 
the deviations {Df^^ „^,Z3^^ „^ }. The exponential decay 
up to round-off error is a strong indication of the regular- 
ity of the solutions, i.e., for the regularity criterion (351. 



An occurrence of logarithmic terms would have become 
apparent by a mere power law decay. 

Note that the perturbation for calculation No. 1 is 
large, as r-u = 2M/a-u « 2.14 rh. The corresponding 
potentials {(j),V^ ,V'^ ^V"^) are displayed in figure [2] In 
contrast, the perturbations are small for the examples 
Nos. 2 and 3. In general, we can assure that our inner 
boundary is an apparent horizon if we choose the per- 
turbation to be small enough. However, for large pertur- 
bations the outermost marginally trapped surface may 
be located further out. To clarify this point, we plan 
to apply a numerical apparent horizon finder in a future 
project. In this paper we stress the fact that we are able 
to construct arbitrarily small perturbations, even in a 
tiny vicinity of extreme Kerr black holes (see calculation 
No. 3). Note that we are not limited by some critical 
value jcrit below the extreme cascj^ 

A careful subsequent analysis of our solutions shows 
that the data constructed indeed correspond to perturbed 
Kerr black holes. In order to verify this point, we com- 
pute at the inner boundary mass and angular multipole 
moments Af„ and J„, which were introduced for isolated 
horizons by Ashtekar et al. [25]. We compare the mo- 
ments with a reference Kerr solution possessing the same 
values for Mg and Ji and find that all higher moments 
are different from the corresponding Kerr expressions. In 
this manner we can conclude that our solutions do not 



In the case of conformally flat initial data, such a critical value 
jcrit has been observed, see e.g., [28] where jcrit ~ 0.97. 
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FIG. 1. For the examples listed in (jlj, the maximal deviations 
Dt,,-a^ (red, dashed), (blue, dot-dashed), -DJ;„,„^, 

(green, dotted) and Dn^^n (black, solid) are displayed, see 



with t 



eqns. (50 51 1. The numerical solutions have been calculated 
;he resolutions (no-,n^) = (no-, \na 



represent the unperturbed Kerr solution in some new co- 
ordinates but do indeed describe perturbed Kerr data. 
An extensive description of this test is given in appendix 
IVIIICI 



VII. DISCUSSION 

In this paper we have constructed single black hole ini- 
tial data on compactified hyperboloidal slices on which 
the mean extrinsic curvature K approaches a constant 
asymptotically (ACMC-condition, K — Kq + Kj^a'^). 
In particular, we were able to identify ACMC-slices for 
the Kerr solution and generated perturbed Kerr data 
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FIG. 2. The potentials {(f), V^} for strongly perturbed Kerr 
initial data are displayed, see example No. 1 in The rota- 
tion parameter is j = 0.5, and the inner boundary is located 
at a-H = 0.5 which corresponds to rn ~ 2.14 rh. 



from these. We find rapid exponential convergence of 
our pseudo-spectral solutions indicating the regularity of 
the potentials. Therefore, we believe that the stronger 
global CMC-condition, i.e., K4 = 0, which is adopted by 
several authors, can be replaced in favour of the local 
ACMC-condition without losing regularity. Note that 
the freedom provided by the ACMC-condition can be 
used to optimize the numerical representation. In this 
regard, a comparison of CMC- and ACMC-slices of the 
Schwarzschild metric is given in Appendix | VIII 5] show- 
ing that the CMC-condition needs more coefficients for 
an accurate spectral approximation. 

The numerical solutions calculated in this work are 
meant to serve as initial data for a time-evolution code 
that solves the dynamical Einstein equations. For the de- 
velopment of such a code, we plan to realize the concept 
suggested by Zenginoglu [5D], i.e., we will treat the hy- 
perboloidal evolution problem of the Einstein equations 
in terms of a conformally rescaled metric. Note that in 
Zenginoglu's approach, the wave gauge is being utilized 
so that the Einstein equations become a system of wave 
equations. By taking as a starting point initial data cor- 
responding to slightly perturbed Kerr black holes, we can 
explore the new code in the weakly non-linear regime 
where the equations should be similar to corresponding 
linearised ones. Numerical methods to realize this project 
are under construction. 
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VIII. APPENDIX 

A. The coordinate transformation leading to 
ACMC-slices 

In section |TT] we described the coordinate transforma- 
tion ([4] |5| which leads, for the Kerr solution, from Kerr 
coordinates to compactified hyperboloidal slices on which 
the ACMC-condition ^ is vahd. The equations ^ [s]) 
contain an auxiliary function A = A{a, jj) which was de- 
termined through an appropriate diagonal Pade expan- 
sion. Here we provide the explicit form of this function: 



with 

Z = 24ci (ciC3 - C2) - {aic\ + 2Acl - 48C1C2C3) a 
N = 2A (ciCg - cfj + (24c2C3 - 0401 )cr 
+ {aiC2 ~ 24c^) 0-2. 
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The coefficients = Ci(fi ; j) are given by: 



ci = ai 
1 

C2 = 7^ i 



16 

y - 4) 



C3 - - ^ (8 + 12ai + 6a? - j") + ^(1 + ai)j V 
6 ^ '8 

■4 4 

256 ' 

To fix the free parameters oi and 04 as functions of j, 
we analyse the decay of the Chebyshev coefficients c^"-* 
of the conformal lapse d (see section |ll] and footnote [s] 
therein). We fit their logarithm linearly over the index, 



Fit 



log(M|cf)|) 



A + /2 • k. 



and perform a numerical minimization of /2 with respect 
to ai and 04 using a simplex algorithm. In this way we 
obtain the parameters corresponding to the most rapid 
decay. We observe that the results obtained in this man- 
ner are well represented by a simple quadratic fit, 



a-i 



04 = 4! 



113 9 
Too ^ 500 
161 3 



+ 40^' 
5 



10 ^ 4'-^' 2-^ 



Note the occurrence of terms containing \j\ which is a 
result of the specific quadratic fit we have chosen. In a 
more fundamental treatment one would try to avoid such 
terms, as the Kerr solution should be described by metric 
coefficients that are smooth with respect to j. However, 
for the purpose of this paper it is convenient to utilize 
the above choice, as it realizes with little effort ACMC- 
slices possessing smooth functions d with rapidly decay- 
ing Chebyshev coefficients. 



B. Comparison of ACMC- with CMC-slices 

In this section we perform a numerical comparison of 
our ACMC-slices with CMC-slices in the case of the 
spherically symmetric Schwarzschild black hole. The 
comparison consists again of analyzing the decay of the 
Chebyshev coefficients cj^^ of the conformal lapse d in 
the two situations. 

For the Schwarzschild metric we have j — 0, and 
the corresponding auxiliary function A, appearing in the 
coordinate transformation leading to ACMC-slices, be- 
comes independent of /i, 

^(''')acmc — 

3(44789179000-2 + 1950390283a) 
5178027300 + 14182095000cr + 28592255449cr2 ' 

A plot of this function is provided in figure |3] Also, 
the corresponding mean curvature -f^ACMc is displayed 



therein. Note that the asymptotic value of the mean 
curvature is [see equations (17 1 and (20)] 



AM^l + ai 2M V 71 



0.514M" 



With the conformal factor Q. — <j/AM, the associated 
conformal lapse turns out to be 



M^{a - l)(72yl'(CT)2 + (1 - 2a2) A'{a) + a + 1 

In figure [I] the decay of the coefficients Mcf ^ is 
displayed. 

We now turn our attention to CMC-slices. In the 
Schwarzschild case, the line element ([T]) reduces to the 
Eddington-Finkelstein form, describing ingoing null rays. 
As described in [5Tk33j . the coordinate transforma- 
tion between the coordinates (V, r, Lp) and coordinates 
{T,r,'d,(p), in which the slices T = constant satisfy the 
CMC-condition, is given by 



T = V 



rdr 

r - 2M 



radr 
(r-2Af)/' 



where 



/ 



1 



2M 

r 



Kr 
"3~ 



(52) 



and K is the constant mean curvature. The parameter C 
represents an additional one-parameter degree of freedom 
in the choice of spherically symmetric CMC-slices in the 
Schwarzschild metric. For slices extending up to we 
have 

8 , 
C > -KM^, 
3 

hence a\r=2M = —f\r=2M < 0. From the relation be- 
tween T and V, we can read off the transformation 
leading from Eddington-Finkelstein-slices to compacti- 
fied CMC-slices. Putting T = 4Mr, we obtain 

2M 



V = AM[t+-- log((T) + Acmc(ct) 

(T 



with the function ^cmc given by: 
f + a 



2a2(a-l)/'^^ + ^°g('^)^a 



Note that despite the apparently singular form, Acmc is 
analytic on the entire interval a £ [0, 1]. 

Using again the conformal factor — a /AM, the as- 
sociated conformal lapse reads 

acMc = ^ v/l44(l - a)a^ + (SKM - 'iCM-^a^y. 
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Now, utilizing a method similar to the one described 
in section VIII A| we find that strong decay rates of the 
corresponding coefficients c^,"^"'^^ are realised if we put 
the free parameters as 



K « 0.32M- 



C « 2.88M^ 



In fi gure [3] the quantities {^acmc: ^cMc} as well as 
the corresponding mean curvatures {i^ACMCj ^cmc}i ob- 
tained for these settings, are plotted. 
Figure |4] shows the convergence of the Chebyshev coef- 



A(a) 

1.5 : 



1 - 




120 



0.5 





FIG. 3. The quantities A = A{a) and MK = MK{a) for 
CMC-slices (dashed) and ACMC-slices (solid) case are plotted 
in the Schwarzschild case j = 0. The slices possess asymptotic 
mean curvature values at J^^ of Ko ~ 0.514Af~^ for the 
ACMC-slicing and Kg « 0.32M-1 for the CMC-slicing. The 
additional parameter C in the family of CMC-slices has been 
set to CA'/"^ ^ 2.88. 

ficients {M cf M cf '^^'^^l. While the coefficients 
|^^^5acmc)| fQj. ^j^g ACMC-slice reach 10-l^ i.e. the or- 
der of machine precision (~ 16 digits), at fc = 32, the 



coefficient iMc' 



(acMc) I 
■32 I 



for the CMC-slice is only about 



10~^^, i.e. larger by 4 orders of magnitude. We thus con- 
clude that the replacement of the global CMC-condition 
in favour of the local ACMC-condition provides advan- 
tages in a numerical treatment. 



C. Mass and Angular Multipole Moments 

Having obtained the numerical solution of the con- 
straint equations, we perform a source multipolar de- 



FIG. 4. Chebyshev coefficients of the conformal lapse a for 
the Schwarzschild metric in CMC-coordinates (dashed) and 
ACMC-coordinates (solid). For the ACMC-slice, the order of 
machine precision, that is 10~^®, is reached at k = 32, while 
for the CMC-slice the corresponding coefficient is of the order 
10-12. 



composition of the marginally outer-trapped surface H 
characterizing the inner boundary of our numerical do- 
main. Although we focus here on small perturbations 
around the Kerr solution, we note that the expansion 
scalar of outgoing null rays on surfaces inside the numer- 
ical domain (cr < a-u) is positive even for the stronger 
perturbation depicted in Fig. ([2]). This result is signifi- 
cant evidence that, also in these cases, we are still dealing 
with apparent horizons, i.e. outermost marginally outer- 
trapped surfaces. 

Introduced in (29j , the multipole moments of isolated 
horizons are an invariant measure on "H that allow us to 
tell if the new initial data represents, indeed, a deformed 
black hole or if it is rather just a new slice of the Kerr 
space-time (with possibly new parameters M and j) |34j . 
In fact, as was shown in [25], the geometry of a given iso- 
lated horizon as well as its vicinity are completely charac- 
terized by multipole modes divided into two types: mass 
modes Mn and angular modes J„ . In the case of dynam- 
ical horizons, they were first employed in j35] , where the 
authors, among other results, used them to verify that 
during the numerical evolution of a space-time, the black 
hole settles down to a Kerr solution. 

In order to define the horizon's angular momentum in a 
gauge-independent manner, the calculation of such multi- 
poles relies on the existence of a preferred divergence-free 
vector field on "H, which here is taken to be the axial 
Killing vector 9^ . Another important issue for their con- 
struction is the introduction of a new coordinate system 
{/i, (f} on H, in such a way that the Legendre polynomi- 
als Pn(A) hold the correct orthonormality properties. As 
described in [2niinS]i S [0,27r) is the affine-parameter 
along which coincides, in our case, with the original 
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coordinate ip. Apart from that, /t = cos9 G [—1,1] is 
defined by 



(53) 



with ^^^Efea the area-element form on Ti. and Ru the areal 
radius defined by the area of %, i.e., A-^ — §^dA — 



The integration constant in ( 53 1 is determined 
by requiring fid A — 0. 

The angular and mass multipoles are then defined as 



Jn 



R 



H 



M„ = M 



Stt j-h 

TDTl 



H 



Sir 



KMP'n{Ki^))dA 
^^^R{fi)Pn{Kf^))dA. 



(54) 
(55) 



n 



In Eq. (541, K^, 



and the derivative of P„ 
should be taken wi th r espect to its argument fi. Apart 
in Eq. (IsHl, 



from that, 



^^^i? is the Ricci scalar of the 
two-surface Ti. and the parameter M-u is given by M-u — 
^JR^/4, + Jl/ R^, which corresponds to black-hole mass 
in case of the Kerr space-time. 

Using the Gauss-Bonnet theorem, it is not hard to con- 
vince oneself that Mq = M-^. One can also show that 
Ml = [ini IMj- Regarding the angular modes, Jq is 
clearly zero and Ji is exactly the definition of the Komar 
angular momentum on the horizon. Finally, due to the 
symmetry across the equatorial plane presented by our 
initial data, all the odd Mn as well as the even J„ vanish. 

Given the numerical initial data, calculated as de- 
scribed in the previous sections, we evaluate the mul- 
tipole modes ( 54 1 and ( 55 1 with respect to the coordi- 



nate sphere cr^. If our solution were just a new slice of 
a Kerr space-time, the corresponding mass and angular 
momentum parameters of the metric would be given, re- 
spectively, by Mq and Ji . In this way, the comparison of 
higher modes allows us to distinguish a perturbed black 
hole from the Kerr solution. 

It is worth to mention that both multipoles modes, as 
defined in eqns. (54 1 and (55 1, have physical dimensions 
of [Mass]"'"''^. However, as we do not want to attach any 
strong physical meaning to the multipoles, but rather 
just use them as an unambiguous tool to discriminate 
our initial data form the Kerr solution, we present our 
results in terms of the normalized modes Jn/MQ^^ and 
M„/Mo"+i. 

Fig. [5] depicts results for small perturbations around 
the Kerr solution. In particular, the top, middle and 
bottom panel have the parameter j = 0.5, j = 0.9 and 
J = 0.999, respectively. The corresponding horizon were 
deformed to a-u — 1.07, cr-u — 1.4 and a-^ — 1.9145 
(their original values are given in Table g. We find that 
the multipole modes confirm our expectation that the 
new initial data represent deformed and rotating black 
holes. In fact, the inset within each panel shows that 
the multipole modes of our perturbed initial data do not 
agree with the ones corresponding to the Kerr solution 
with the same parameters Mq and Ji . 
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FIG. 5. The top, middle and bottom panel show, respectively, 
the multipole modes for solutions obtained with j = 0.5, 
j = 0.9, and j = 0.999 and inner boundary set, respectively, 
to a-H = 1.07, (Tfi = 1.4, a-H = 1.9145 (original values are 
given m Table The non-agreement between them confirms 
that our initial data do not correspond to a different slice 
of the Kerr solution. Crosses and dots correspond, respec- 
tively, to the regularized mass \Mn\ /Mg+'^ and the angular 
\Jn\ /Mq'^^ modes, while the inset shows the difference be- 
tween their values for the perturbed black hole and the corre- 
sponding Kerr solution with the same Mq and Ji parameters. 
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